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Abstract 

A semi-empirical model of the infrared (IR) spectrum of the Cgo molecule is 
proposed. The weak IR-active modes seen experimentally in a Ceo crystalline 
sample are argued to be combination modes caused by anharmonicity. The 
origin of these 2-mode excitations can be either mechanical (anharmonic in- 
teratomic forces) or electrical (nonlinear dipole-moment expansion in normal 
modes coordinates). It is shown that the electrical anharmonicity model ex- 
hibits basic features of the experimental spectrum while nonlinear dynamics 
would lead to a qualitatively different overall picture. 

PACS numbers: 33.20.Ea, 33.20.Tp, 33.70.-w,78.30.-j 
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I. INTRODUCTION 



There has been a great deal of progress in our understanding of the chemistry and 
physical properties of fullerenes. Discovery of superconductivity in alkali-metal doped Cqq 
[ l| has ignited discussions on possible mechanisms of this phenomenon One class of 

models stresses the coupling between electrons and intra- molecular phonons ||. Raman 
and infrared (IR) spectroscopy have probed the vibrational properties of Cgo compounds 
P^TDf and many theoretical models have tried to explain properties of the 46 distinct modes 
predicted by group theory. 

The icosahedral (Ih) symmetry of Ceo allows four distinct IR active modes (Ti u ) and ten 
Raman active modes (2A g 8H g ) in harmonic approximation. It is customary to denote the 
IR modes at frequencies 528, 577, 1183, and 1429 cm -1 , as T lu (i), i=l,2,3,4, respectively. 32 
optically inactive (silent) modes are 1A U , 3T lg , 4T 2g , 5T 2u , 6G g , 6G U , and 7H U . Higher order 
peaks are seen experimentally by increasing the optical depth of a sample. In principle there 
are 380 second-order combination modes IR allowed by the Ih symmetry ||. Second-order 
overtones are IR forbidden. 

Several authors reported observation of weak modes in Raman [|i~0| , |i~Tf and IR |7|-|9|JT2]| 
spectroscopy. Wang et al. 0, Martin et al. ||, and Kamaras et al. analyzed the weakly- 
active features in conjunction with Raman |TtJ and neutron measurements [|T3] to extract the 
32 fundamental frequencies of the silent modes. The frequencies differ significantly among 
the authors leaving the question of the assignment of fundamentals open. 

Possible mechanisms of activating the weak modes include 13 C isotopic impurities, crystal 
environment effects and anharmonicity. Impurities, dislocations and electric field gradients 
at surface boundaries can be excluded due to their sample dependence. An experimental 
and theoretical vibrational study of 13 C-enriched crystals excluded the isotopic symmetry 
breaking as a potential candidate |14| . A few of the weak modes are thought to be activated 
due to the fee crystal field effect. The crystal field reduces the Ih symmetry of Ceo and 
activates silent odd-parity modes. Above 260K the C 60 molecules freely rotate and the 
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time averaged crystal field perturbation is zero. This effect of 'motional diminishing' of 
silent modes have been experimentally observed and theoretically studied by Mihaly and 



Martin An experimental study of pressure dependence of these modes would help to 
substantiate this mechanism. 

The goal of the present paper is to identify, qualitatively, the mechanism of activation 
of the higher-order vibrations; detailed assignment to normal modes remains a task for the 
future. Basic formalism of anharmonic effects on IR activity is given in |T^-^TJ. There 
are two ways in which anharmonicity can display itself in an optical spectrum. Either it is 
driven by anharmonic interatomic forces (mechanical anharmonicity) or by an anharmonic 
coupling of a photon field to two or more phonons (electrical anharmonicity). Although the 
two mechanisms are not independent, each has its own characteristic absorption intensity 
pattern. When compared with an experimental spectrum one can decide which of the two 
kinds of anharmonicity prevails in the IR spectrum of Cgo- Although the spectrum may 
contain cross-contributions from both phenomena, here they are treated separately. 

Several models have been used to calculate absorption intensities in harmonic approxi- 
mation. Tight-binding models |22||2B| are in complete disagreement with the experimental 
results. The bond-charge model [23] fits very well with frequency positions of fundamentals 
but the IR intensity pattern disagrees with basic trends in the observed spectrum. The 
same is true for a Hubbard type model stressing electronic correlation effects |[25|| . Relative 



intensities are best reproduced by the LDA approximation |23|,|26j . Due to its computational 
complexity the LDA scheme is not convenient for computing second order intensities. We 
therefore propose a semi-empirical model which is satisfactory for a qualitative comparison 
with experiment. Figure 1 summarizes the performance of these models in calculating the 
absorption intensities. 

Some characteristics of the experimental IR spectrum || are shown in Fig. 2. Combi- 
nation (difference) modes are higher-order modes with frequency u equal to Ui ± Uj, the 
sum (difference) of fundamental frequencies tUj. Their intensities are temperature dependent 
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according to (nj + |) ± (nj + |), where Hi is the Bose factor, + | = | cot^^l^), with a 
temperature T and the Boltzmann constant Following features can be observed in the 
spectra: (i) besides four first-order peaks there are more than 180 weak absorptions, (ii) no 
difference peaks are resolved (i.e. no temperature dependence of intensities except a trivial 
improvement in the frequency resolution at lower temperatures), (iii) most of the spectral 
weight is in the high-frequency regime (1000 - 3000 cm -1 ), and (iv) weak modes around four 
first-order bands are not enhanced through a resonance effect. 

This paper treats the frequency positions and absorption intensities independently. Nor- 
mal modes and frequencies are calculated using a simple force-constant model proposed by 
Weeks [p7| . This model fits IR data reasonably well but is not expected to give especially 
realistic eigenfrequencies for the silent modes. The dipole moment which arises due to the 
electron-phonon coupling determines the absorption intensities |28f . Only second-order com- 
bination and difference modes are considered in the paper. Sec. || deals with the mechanical 
anharmonicity problem with the Morse function used for the interatomic bond-stretching 
potential f27fl . A linear relation between the dipole moment and ionic coordinates is pro- 
posed in this section. The relation contains parameters fittable to the relative harmonic 
absorption intensities. Second-order modes are computed using a perturbation method ig- 
noring possible resonances. However, the intensity pattern of the second-order modes fails 
to reproduce experimental features. An electrical anharmonicity model is therefore intro- 
duced in Sec. [TTI|. Normal frequencies and normal modes are again taken to be those of 
the Weeks' model. A semi-empirical model for an electronic configuration on a distorted 
Ceo is presented which allows the electronic coordinates to depend in a nonlinear fashion 
on positions of ions. This gives rise to an intensity pattern very similar to the experimental 
one. Finally, conclusions are drawn in Sec. |TV|. 
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II. MECHANICAL ANHARMONICITY MODEL 



Considering the C 60 molecule as a system of oscillating ions with electrons moving adi- 
abatically in their field, the ionic dynamics is governed by the potential: 

46 Si y 46 9i,9j,9k 

V — 2 1212 mU i Qiq + 7T 12 12 CigJr^sQigQjrQks- (1) 
i=l q=l i,j,k=l q,r,s=l 

Here Q iq is the g-th normal mode coordinate belonging to the frequency u>i, i—l,..,46, 
q=l,..,gi, and g; t is the degeneracy of the i-th band. Higher-order terms are neglected. 
The anharmonicity coefficients Ci q j r ^ s are given by 

C ^ (2) 

q,J ' dQi q dQj r dQks 

Light couples to the system via the term 

V 1 = -/i(Q) • E, (3) 

where E is the externally applied macroscopic electric field and ji stands for the dipole 
moment of the system. The latter is generally a nonlinear function of normal coordinates 

46 9i ^ 46 9j,9k 

^ = 1212 M igQig + 2 12 12 M jrM Q jr Q ks (4) 

i=l q=l j,k=l r,s=l 

Again, higher-order terms are not included and the folowing formulas determine the expan- 
sion parameters M ig and M jrj fc s : 

and 

J ' dQj r dQks ^ ^ 

The vectors Mj g are nonzero only when the Qi q mode is IR allowed. 

Anharmonic dynamics {Ci q j r ^ s ^ 0) and a linear coupling of light to phonons (M^^ = 
0) characterize the mechanical anharmonicity (MA) phenomenon. 
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Several force-constant models for Ceo have been presented [^7],^|3(]]. To calculate normal 



coordinates and the anharmonicity coefficients Ci q j rj k s , I use the model suggested by Weeks 
p7| which is a refined model of Weeks and Harter . This model contains two parameters 
which were fitted to selected IR and Raman frequencies. Ionic dynamics is governed by two 
types of interactions: (i) the Morse potential producing anharmonic terms 

90 

V m = Y,D{1 - explain - r eq )}} 2 , (7) 
i=i 

controls bond-stretching. Here D, a, r eq , and rj are, respectively, the dissociation energy, 
Morse anharmonicity, equilibrium and instantaneous length of the ith bond. Summation 
runs over all bonds. The dissociation energy is estimated as the average of the dissociation 
energies of a single and a double C2 bond, D = 5.0eV, the equilibrium length is taken to be 
1.4A and the parameter a was fitted to the value 1.6A \ (ii). The bond-bending harmonic 
potential is given by 

Vb = J2 7 l^eq-0j) 2 , (8) 

j 

where the summation is over the 60 pentagonal angles with the equilibrium angle of |7r and 
the 120 hexagonal angles with the equilibrium angle of |7r. The potential does not distinguish 
between hexagonal and pentagonal angles and the best fit yields 77 = 12.48eV/rad. 

The bond-stretching potential in the harmonic approximation together with the bond- 
bending potential give normal coordinates and frequencies. The coefficients Ci q ^ r ^ s come 
from the expansion of the Morse function to the third order in ionic distortions from equi- 
librium and from the transformation of the Cartesian coordinates to the normal mode ones 
computed numerically. Qualitative behavior of the normal modes of the model (with the 
bond-stretching potential in the harmonic approximation) is discussed in the original pa- 



pers [2^,3J]. It is enough to note that lower frequency normal modes exhibit mostly radial 
distortions while the motion of higher-frequency ones is tangential. 

IR intensity of a given mode is proportional to the square of a dipole moment associated 
with the mode. If ionic charges of the same value were put on the vertices of Ceo, the 
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resulting dipole moment would be zero due to the center-of-mass conservation. The dipole 
activity is therefore caused by changes in the electronic configuration. Carbon valence 
electrons fall into two classes. The first class consists of cr-electrons positioned with the 
highest probability in the middle of bonds. These electrons have fixed charges and do not 
contribute to the dipole moment (due to the center-of-mass conservation). In the following 
the notion of a bond charge will include also a contribution from ions in some effective way. 
The sign of such an effective bond charge will not be important, it can be either positive 
or negative. Allowing the bond charges to acquire a charge with dependence on the bond 
lengths or by some other mechanism leads to a spectrum where the Ti u (2) mode is hardly 



visible instead of having the second largest activity p^ , |32| . The second class consists of 
7r-electrons which create a dipole moment in the following way. Consider these 7r-electrons 
to be vertex electrons moving in the field of their parent ions. Let these electrons interact 
further only with the three nearest ions. The positions of the 7r-electrons are modeled in the 
following way. Denote the radius-vector of the i-th electron measured from the vertex i 
with the position Rj and R^, j = 1, 2, 3, the nearest ions positions, respectively, seen from 
the center of Cgo- The direction of is taken to be the direction of the normal vector rij to 
the plane given by three nearest ions with a rescaled position of the one making the double 
bond with the vertex. This condition, 

n, • (R? - R«) = n, ■ (R» - Cl R?) = 0, (9) 

introduces a fitting parameter ci, effectively measuring the ratio of the double- and single- 
bond charge (here the bond Rj — R^ is the double-one). Single bonds are bonds connecting 
a hexagon with a pentagon and double-bonds are connecting two hexagons. When there is 
more charge on the double-bond than on the single-one, the parameter c\ is greater than 
unity. If the bond charge is negative, the direction is out of the bucky-sphere and if it is 
positive, the direction is inwards. 

Consider the distance d of the vertex ion to the plane given by its three nearest ionic 
neighbors (with the double-bond neighbor rescaled as explained above). Denote as d eq the 



distance for the equilibrium configuration. Let, for a moment, the effective bond charge be 
negative. If a distortion of the ionic positions occurs such that d > d eq the vertex electron 
will be pushed 'out' of the Ceo sphere and vice versa. If the net bond-charge is positive, 
the situation is inverse. This phenomenology reflects a Coulomb repulsion (attraction) of 
the vertex electron by (to) adjacent bonds. When these bonds move closer together the 
vertex electronic cloud is deformed such that the mean electronic position will be as far 
(close) as possible from (to) the bonds. The effective rate of the deformation will be the 
second free parameter c 2 (the same for each vertex due to symmetry). The relation between 
the electronic position and the distance between the vertex ion and the plane given by its 
nearest neighbors can then be expressed as follows : 

r t = {1 + c 2 [d i (c l ) - d e9 (ci)]}nj(ci), (10) 

where the dependence on the parameter c\ is indicated. The dipole moment is then clearly 

60 

y. = ^[l + c 2 {di - 4<?)]ni. (11) 

i=i 

The normalization in both formulas is not important for calculating relative values. The 
distances di depend for small distortions linearly on normal coordinates, so only the linear 
term is kept here because the mechanical anharmonicity couples this linear displacement to 
two normal modes. 

There are two natural parameters in this model, c\ and c 2 . In the harmonic approxima- 



tion the intensity of the j-th mode is [28 



1$ = £ M %- (12) 

3=1 

Experimentally obtained relative intensities are 1., 0.48, 0.45, 0.378 for the modes T lu (l), 
T lu (2), Ti u (3), T lu (4), respectively P^j . The best fit to these intensities yields the values 
Ci = 1.59 and c 2 = 0.67A . The IR spectrum obtained with the fit (all peaks in this and 
following figures have the Lorentzian widths taken to be uniformly 2cm -1 ) along with an 
experimental one is shown in Fig. 3. Agreement with experiment is very good. 
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For the frequencies that are not in the immediate neighborhood of the frequencies of 
the four IR allowed fundamentals, the following formulas were obtained in Ref. 18 for the 
second-order intensities of combination and difference modes: 

Ct,= sfer^U + nfc + m) (13) 



9k,9l , 

£ ( I 

r,s=l jElR 



9 *4 l ( v (MjICj,^,,,) n 2 



and 



v (M. 3 \C jMs ) A2 
r,s=l fc; 

respectively. The summation in brackets is over four IR active bands and the inner-product 
notation stands for the sum over a degenerate set : 

(M,0...) )►>,/:,,...• (15) 

When the frequency of a combination (difference) mode is near the frequency of an IR 
allowed mode (the Fermi resonance effect), a perturbation leads to a mixing of the two modes 
and spreads out their frequencies (see Ref. [|HJ). The second-order modes are enhanced 
conserving the original spectral weight so the integrated absorption intensity of the band 
is unchanged by the anharmonic perturbation. If the spectral resolution is not enough 
to resolve the two modes the resulting picture is similar to the original one without a 
perturbation. The Fermi resonance effect has not been observed in Cgo- 



Figure 4 shows the results of the numerical calculations based on the Eqs. |l3^|14] . Some 
trends in the spectrum are clear already from the equations. First of all the second-order 
intensities are relatively weak compared to the experimental spectrum in Fig. 2 (the ex- 
perimental picture here is somewhat misleading due to the saturation of first-order peaks). 
Most intense modes have frequencies close to the four IR bands, leaving high-frequency com- 
bination modes practically invisible. Moreover there are relatively intense difference modes 
(identified by their strong temperature dependence) in the lower part of the spectrum. These 
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features are in contradiction to experiment thus excluding mechanical anharmonicity as the 
mechanism for activation of the combination modes seen in experiment. In matching the 
combination modes to experimental data, authors in Ref. 8 did not find any evidence for 
a significant deviation of the frequencies of these modes from the values of u>i + uij. This 
supports the above conclusion that mechanical anharmonicity is not producing significant 
effects in the Cqo IR spectrum, since the relative frequency shift as a consequence of me- 
chanical anharmonicity only is of the same order of magnitude as the relative intensities of 
the second-order modes. 



The electrical anharmonicity (EA) is a less studied phenomenon of molecular physics 
than the mechanical one. It is based on the fact that the dipole moment is generally a 
nonlinear function of normal modes. In view of Eqns. [I] and |], electrical anharmonicity 



Selection rules for the second-order modes are reflected in the elements of the matrix ~M.j r ^ s , 
and are the same as in the case of the mechanical anharmonicity. Since the ionic dipole 
moment is linear in ionic positions it is clear that the nonlinear contribution stems from a 
nonlinear response of electronic positions to a change in ionic configuration. A harmonic 
treatment now suffices for the ionic displacements; the Weeks model of Sec. II. is used. 

The nonlinear electronic response is modeled in the following way. The notation is the 
same as in the previous section. Consider again a 7r-electron in the field of its parent ion 
and adjacent bond-charges. The interaction with its nearest neighbor ions is governed by 
the Coulomb potential 



III. ELECTRICAL ANHARMONICITY MODEL 



arises from the second term in Eq. f|, while the ionic dynamics is harmonic (C, 



;ks— 0)- 




(16) 



and similarly the interaction with adjacent bond-electrons is given by 




(17) 
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Summations are over the three nearest ions and Rj is the position of the vertex ion. Note 
that while R's are measured from the mass center of Cgo, r, is measured from the position 
of the ith vertex ion (Rj). The strengths of the interactions are measured by some effective 
charges K{ and Kb e for neighbor ions and adjacent bond-electrons respectively. Only the ratio 
— is a relevant fitting parameter. The motion of the 7r-electron in the field of its vertex ion 
is simplified by restricting it to a sphere around the ion with a radius R which will be the 
second fitting parameter: 

r, = Rn t . (18) 

This gives a simple two-dimensional minimization scheme: for each vertex and a pair of 
fitting parameters (R, fti/fi^e) find a unit vector n, such that the function 

K-iO^+Ve-feeK) (19) 

is minimal. The electrical dipole moment is then computed and resulting first-order inten- 
sities (Eq. |12|) compared with corresponding experimental values. The best fit corresponds 
to values of R = 0.06A and ki/k^ = 4.80. For some range of the parameters there are two 
electron positions for which the potential in Eq. |H] has a local minimum. In such cases the 
global one was considered. The best fit lies in the region with one minimum. It is obvious 
that the best fits have no physical justification. To support the model I did simulations with 
different, more physical values of the free parameters obtaining the same qualitative picture 
as will be shown later. It is also appropriate to remark that a feedback from the adiabatic 
changes in electronic positions to ionic motion is implicitly considered in the harmonic level 
in the force-constant model. 

For the IR absorption the changes of the minima positions with ionic distortions are 
relevant. Numerical differentiation was used to obtain the dipole-moment matrices and 
ISAjr^g from Eqns. |5] and |6]. An important feature of the model is that the 7r-electronic 
positions are more sensitive to tangential distortions than to radial ones. 

Second-order absorption intensities of combination and difference modes have now simple 
forms [ |l~8|l : 
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TEA _ ^ U k + U)i f , 9 M l 2 



"a™ ^ > 1 + "' + "')£ ML - (20) 



r,s=l 



Figure 5 shows the spectrum obtained from equations |2CJ and [21]. The following features 
can be extracted. Overall intensity of the weak modes is higher (in relative sense) than 
in the case of the mechanical anharmonicity. Spectral weight is shifted towards higher 
frequencies. This is a consequence of high sensitivity of electronic positions to tangential 
distortions which are characteristic for higher-frequency modes. The sensitivity of electrons 
to the tangential ionic motion is also the reason why difference peaks have relatively very 
small intensity (the difference peaks are most intense in the region of 600 - 1000 cm -1 , 
however the intensities are much smaller than those of combination modes in the region 
1000 - 3500 cm -1 ). There is obviously no resonance effect since the two terms in the Eq. |] 
are independent. The frequency distribution in the Weeks model differs from that in Ceo so 
a closer comparison with experiment is not possible. One consequence is that in Figure 5 
weak features up to 4000cm -1 are visible, while experimentally weak peaks above 3500cm -1 
have not been resolved. This difference in the frequency distribution may be a part of the 
reason why there is so little activity in the region 600 — 1000cm -1 . Note that almost all of 
the peaks experimentally observed in this region were associated with modes IR forbidden 
in the second-order [|7|,|8| and their appearance must be accounted for by other mechanisms. 



IV. CONCLUSION 

Mechanical and electrical anharmonicity provide possible mechanisms for activating weak 
modes resolved in IR spectra of Ceo thin films and single crystals. I have proposed simple 
semi-empirical models of the phonomena. The main features of the models are: (i) separation 
of ionic dynamics and mechanism of optical activation (the models can be used for any set 
of normal modes), and (ii) emphasis on 7r-electronic system rather than on bond charges. 
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Both models give a spectrum of combination and difference modes which is compared with 
IR measurements. It is found that mechanical anharmonicity exhibits features different than 
those observed. These features can be generally expected from basic formulas (e.g. those 
of Eqs. [L3| and 0) and the model described in Sec. II only helps to visualize them. As a 
by-product the intensities of four first-order IR allowed bands are well reproduced. 

The electrical anharmonicity model introduced in Sec. Ill is based on a nonlinear re- 
sponse of 7r-electronic configuration to ionic distortions. Now the absorption spectrum has 
fewer characteristics given a priori by a theoretical formula and is more model-dependent. 
The main feature the model which leads to quite successful comparison of its spectrum with 
experiment is that electronic positions are much more sensitive to tangential ionic motions 
than to radial ones. 

The separation of mechanical and electrical anharmonicity is posteriorly justified by the 
dominance of the latter. However, the IR activity around four first-order peaks is caused by 
mechanical anharmonicity due to resonance effects, as discussed in Sec. II. There is still a 
region of optical activity (600 — 1000cm -1 ) which this simple model cannot explain. Although 
trial assignments exclude most of the observed peaks in the region as combination modes, 
the question is still an open one and more sophisticated quantum-mechanical treatment can 
yield more authoritative results. 
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FIGURES 

FIG. 1. Comparison of calculated relative absorption intensities of IR allowed T lu (i), 
i=l, 2, 3, and 4, modes with experiment. Intensities of the band Ti u (l) are taken to be 
unity. 

FIG. 2. Ceo single crystal IR transmission spectra at 300K and 77K by M. Martin et. 
al. ||. 

FIG. 3. First-order IR allowed intensities calculated in Sec. II and experimentally 
obtained spectrum (inset) by Hare et.al ||. 

FIG. 4. IR spectra at 300K and 77K computed using the mechanical anharmonicity 
model introduced in Sec. II. Difference modes are easily identified by their strong tempera- 
ture dependence, while combination modes show no such trends. 

FIG. 5. The electronic anharmonicity model (Sec. Ill) produces absorption spectra 
which show similar trends as experimental ones. Difference peaks carry very little spectral 
weight comparing to high-frequency combination ones. 
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□ This work. Sec. II. 
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